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In this paper we investigate a steady accretion within the Einstein-Straus vacuole, in the presence 
of the cosmological constant. The dark energy damps the mass accretion rate and — above certain 
limit — completely stops the steady accretion onto black holes, which in particular is prohibited in 
the inflation era and after (roughly) 10 12 years from Big Bang (assuming the presently known value 
of the cosmological constant). Steady accretion would not exist in the late phases of the Penrose's 
scenario - known as the Weyl curvature hypothesis - of the evolution of the Universe. 
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I. INTRODUCTION 

The classical Bondi accretion model [T] describes, in 
the test gas and stationarity approximations, the spher- 
ical accretion/wind of a barotropic fluid onto a newto- 
nian gravity center. That description has been later ex- 
tended to general-relativistic spacetimes. The steady ac- 
cretion in the Schwarzschildean geometry has been in- 
vestigated by Michel [5] and Shapiro and Teukolsky [3J, 
and in a spherical symmetric spacetime with backreac- 
tion by Malec [3]. A remarkable universal behaviour has 
been discovered for selfgravitating transonic flows, un- 
der suitable boundary conditions: the ratio of respective 
mass accretion rates Mgr/Mn appears independent of 
the fractional mass of the gas and depends only on the 
asymptotic temperature. It is close to 1 in the regime of 
low asymptotic temperatures and can grow by one order 
of magnitude at high temperatures [S]-[B]. Here Mqr and 
Mjv are the general relativistic and the Newtonian mass 
accretion rates, respectively. The stability of steady ac- 
cretion has been established by Mach and his coworkers 

m-m- 

In this paper we consider the accretion onto a black 
hole that is immersed in a cosmological universe. We 
adopt the Einstein-Straus "swiss cheese" model [TU]. A 
spherically symmetric black hole surrounded with accret- 
ing gas sits inside a vacuole that is matched — through 
the D amour-Israel junction conditions [TT], — to an 
external FRLW spacetime. We calculate how the pres- 
ence of a cosmological constant A would affect the most 
important system characteristic — the mass accretion 
rate M of a fluid onto the center. We found that the 
mass accretion rate is damped by the presence of the 
cosmological constant. This effect is stronger for nega- 
tive values than for positive values of the cosmological 
constant and it depends on the ratio of the dark energy 
density to the fluid density. 

In the next two sections we remind the Kottler solu- 
tions and draw the general picture of the steady accre- 
tion, together with suitable equations. Sections 4 and 5 
present three analytic sets of results. These show that the 
matter distribution can be affected by the cosmological 
system. More importantly, we prove that if the cosmo- 



logical constant exceeds certain limit, then the steady 
accretion ceases to exist. Section 6 presents results of 
numerical integrations. It appears that the presence of 
dark energy, either negative (attractive) or positive (re- 
pulsive), damps the mass accretion rate. It appears - in 
accordance with analytic predictions - that the steady 
accretion does not exist above certain values of the cos- 
mological constant. The last section summarizes main re- 
sults and applies them to two cosmological epochs when 
the dark energy is dominant. It is pointed out that our 
results do not agree with the Penrose's " Weyl curvature 
hypothesis" . 



II. ACCRETION IN AN EINSTEIN-STRAUS 
VACUOLE 

We assume relativistic units with G = c = 1. Let a 
FLRW spacetime be filled with dust and dark energy. 
Its line element reads 

ds 2 = -dr 2 + a 2 (r) (j^j-i + r 2 cKlA , (1) 

where the variables have obvious meaning [13] and k — 
0,1 1. We do not need to analyze the FLRW equations, 
for a reason explained below. 

We cut off, following Einstein and Straus [TU], a ball 
from this FLRW spacetime and insert instead a spherical 
symmetric inhomogeneity — a black hole surrounded by 
a spherical cloud of gas — satisfying Einsteins equation. 
It is known that the boundary of the ball must be co- 
moving if the FLRW geometry contains only dust and 
the dark energy represented by A, assuming the absence 
of a boundary layer [T3], [T5] . 

Inside the excision ball itself we can distinguish an in- 
ner accretion region such that the areal radius R does 
not exceed a limiting radius Rqo (that in turn should be 
smaller, or perhaps much smaller, than the areal radius 
of the boundary of the ball), and the asymptotic vacuum 
region outside i?oo- In the asymptotic region the metric 
is given by the Kottler ( Schwarzschild-de Sitter (SdS, 
A > 0) and Schwarzschild-anti de Sitter (SAdS, A < 0) ) 
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spacetime line element 



ds 2 = - 1 - 



2m A 



R 



R z dt 



dR 2 



1 - 



- #i? 2 



-R 2 dn 2 . 

(2) 

The Darmois-Israel gluing conditions [TT], [T^j demand 
that along the comoving boundary the first and second 
fundamental forms are continuous. In the spherically 
symmetric case the radius R has to be continuous and 
in addition: i) the excess mass m should be exactly 
equal to the mass of a dust in the excised ball; ii) the 
boundary should be comoving with the Hubble velocity 
H. Moreover, the stress-energy tensor T^ v should be 
continuous across the junction surface. Balbinot et al. 
have found an explicit solution that describes the corre- 
sponding matching between the Schwarzschild-de Sitter 
and FLRW (dust+A) spacetimes [T7] . Therefore we stop 
discussion on FLRW and SdS (SAdS) solutions at this 
point and focus on the description of the accretion zone. 

We assume that in the accretion region the Einstein 
equations can be approximated by a set of stationary 
equations in time intervals that are much smaller than a 
characteristic time T that is defined in the next section. 
Let Roo be the size of the cloud of gas and the asymp- 
totic speed of sound. Outside R^ the geometry can be 
connected smoothly to the SdS or SAdS spacetime ge- 
ometry by a transient zone of a negligible mass. Let us 
remark that the assumption of the approximate station- 
arity can be checked aposteriori, after finding appropriate 
solutions. Approximately stationary accretion equations 
have been derived in 0]. A quasi-stationary solution, 
of a similar accretion problem in the asymptotically flat 
spacetime with a spherical black black hole, appears to 
be stable under axially symmetric perturbations [5]. 



III. EQUATIONS OF STEADY ACCRETION 

The metric in the accretion zone is spherically sym- 
metric and has the form 



ds 2 = -N 2 dt 2 + Mr 2 + R 2 (d9 2 + sin 2 {0)d<f> 2 



(3) 



where we use comoving coordinates t,r,0 < 9 < 7r, < 
<t> < 2-7T — time, coordinate radius and two angle vari- 
ables, respectively. R denotes the areal radius and N is 
the lapse, a is the radial-radial metric component. The 
radial velocity of gas is given by U = ^ . 

The energy-momentum tensor reads T^ v = 
{p + p)u^u„ + V9iiv with the time-like normalized 



four- velocity £/,,, UJJ^ 



T. A comoving observer 



would measure local mass density p — T^U^U^. Let 
be the unit normal to a coordinate sphere lying in the 
hypersurface t = const and let k be the related mean 
curvature scalar, k = ^V;n l = ~4rd r R. We assume 

the perfect gas equation of state p = (F — 1) poe, where 
e is the specific internal energy and T is a constant. 
Assuming that the accretion is isentropic, one can 



derive the polytropic equation of state p — KpQ, with 
a constant K. The internal energy E — poe and the 
rest p and baryonic po mass densities are related by 
p = Po + E = p a +p/(T- 1). 

The matter-related and geometric quantities satisfy 
Einstein equations and the baryonic mass conservation. 
(But let us point out that for isentropic flows the conser- 
vation of baryonic mass also follows from Einstein equa- 
tions [3].) One can find the mean curvature k from the 
Einstein constraint equations G^o — SttT^q [1] 



where m{R) is the quasilocal mass, 

m(R) = m — An / drr 2 p. 
Jr 



(4) 



(5) 



In the line element ^ we have the comoving time. In 
the polar gauge foliation one has a new time ts(t, r) with 
dt s = dt — NUdR. The field d ts is tangent to the cylinder 
of constant areal radius, dt s R — 0. 
One can show that 



M p = d ts m(R) = AttNUR 2 (p + p) 



(6) 



The mass contained in an annulus (R, Roc) changes if the 
fluxes on the right hand side, one directed outward and 
the other inward, do not cancel. The baryonic current 
density reads j 1 * = PqU 11 . Its continuity equation reads 

V M j" - 0. (7) 
For stationary flows the local baryonic flux 



M = -4ttUR 2 p . 



(8) 



is time-independent at a fixed R. 

We say that the accretion process is stationary (or 
quasi-stationary) if all relevant physically observables, 
that are measured at a fixed areal radius R, remain ap- 
proximately constant during time intervals much smaller 
than the runaway instability time scale T = M/M. 
That means that d ts X = (d t - NUd R )X = for 
X = p ,p,j,U...,M. 

For quasistationary flows we have dnM p = and 
d R M = [4]. Therefore both M p and M are equal, 
modulo a constant factor; they can be identified. The 
quantity M will be called the mass accretion rate. One 
obtains from (fSp an expression 



d R U 2 = 



AU 2 
~R 



2U 2 d R \n{ P0 ) 



(9) 



The speed of sound is defined as a = Jd p p. It is a 
useful and straightforward exercise to express the hydro- 
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dynamic quantities in terms of a: 
T - 1 a 2 



P = Po 



P = Po 



PO — POoo 




(10) 



Here (and below) quantities with the suffix oo do re- 
fer to their asymptotic values (i. e., at Rao). Notice 
that p, po and a 2 show the same monotonicity behaviour: 
9rPo = Cidnp — C-2dRa 2 , where C\^Ci are strictly pos- 
itive functions. 

There are two conservation equations that originate 
from the contracted Bianchi identities V^Tjf = 0. One 
of them can be eliminated, if we choose instead the bary- 
onic mass conservation. The case when v = is the 
relativistic version of the Euler equation, 



N M P+iP + p) JR N = °- 



One can solve Eq. (|Tlj), using Eq. (TlO 



N = C(T-l-a 2 



(11) 



(12) 



The whole system of algebraic equations Q — ( 12 ) closes 
with the imposition of the Einstein equation, G rr — 
87rT rr , which is the only integro-differential equation. It 
can be put in the following form 

d , / 2 \ r - 1 - a 2 



a 2 - u 



k 2 



1 



'm{R) _ 2 
k 2 R\ R 



4irR 2 p - 



AR 2 



(13) 



Finally the line element, in (t, R) coordinates, is given by 



N 



ds 2 = - (N 2 - U 2 ) dt 2 -2 — dtdR 



dR 2 

1^ 



-R 2 dVt 2 . (14) 



We shall study transsonic accretion flows. For them the 
principal object of interest is a sonic point, where vanish 
both the denominator a 2 — and the numerator 171 ffi — 
2U 2 +4irR 2 p-^f 2 - of Eq. (13). The corresponding value 



of the areal radius is called the sonic radius, denoted as 
i?*. A closer inspection of the sonic point shows that 
it is a critical point, with branching pairs of solutions 
describing accretion or wind. In the accretion branch, 
below the sonic point the infall velocity \U\/k is bigger 
than a, while outside the sonic sphere the converse is 
true. This analysis is very much standard, replicating 
the work done in the Newtonian case by Bondi [T] and in 
the general- relativistic case by Malec [3]. 



IV. QUALITATIVE ANALYSIS OF SOLUTIONS: 
SDS BLACK HOLES 



In this section we shall investigate some properties of 
transsonic solutions for positive values of A. The speed of 
sound a (and thus also po and p) decreases as a function 
of R — for accreting solutions — when the cosmological 
constant is absent. We shall show, In Lemma 2, that for 
large values of AR 2 the converse is possible. 

Lemma 1. 

Assume A < . Then the speed of sound, the pres- 
sure and the mass density po are decreasing outside the 



sonic sphere, j^X < for X 
Proof of Lemma 1 

Let us define C(R) = 
and let C, 



LEI 



,p and po- 



2U 2 (R) 



AR 2 



C(Rqo). One easily obtains from this def- 



inition, using equation Q in order to eliminate OrU 2 , 
that 



dR 



C(R) 



-AC 
R 



AirR 



6p — 



A 
2^ 



3m o ~rhPo . „o d 
W +4U ^ + ^ R dR P - 



(15) 



The condition A < |j£ 



guarantees that Coo > and — 
using the argument of continuity — the positivity of the 
function C(R) in an open interval (i? s ,i?oo). We show 
that R s = R*. 

Notice that while at the sonic point C(R Sf ) = 0, the 
function 8rC must be strictly positive. Indeed, assume 
the opposite; then from continuity C would have to be 
negative just above R* . Then from Eq. (151 would follow 



that both terms with deri vati ves must be negative, th at is 
— see the remark below ( 10 1 — Orcl 2 < 0. But Eq. ( 13 ) 
gives Bro 2 > 0, if C < and R> R*. This contradiction 
proves our statement and implies the positivity of C in 
an interval (i?*,i?f,) . 

Now, assume that C changes sign at that is 
C{R b ) = and -4C(i?) 



) R =R b < 0- But if C ( R b) = 
then — from Eq. ( 13 ) — the derivative BrO 2 would have 
to vanish. From ( 10 1 would then follow OrPo — 8rp = 0. 
Thus Eq. (fT5l yields at R b 



_d_ 



C = AnR b p (R b ) + 6p (R b ) - 



3m 



(16) 



but this implies j^C > at R b , due to the condition 
A < ySr, This contradiction proves that C(R) is strictly 
positive outside the sonic sphere. 

The positivity of C allows one to conclude — from Eq. 



( 13 1 — that Drci is strictly decreasing outside the sonic 



sphere. This in turn implies the decrease of the pressure 
p and the baryonic mass density po (see the remark under 
Eq. @). 
Lemma 2. 

Assume AR 2 X > ^'g and define R = (^f)*. The 
necessary condition for the existence of solutions is a < 1. 
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If a > 1/3 then R <C Rqo and the speed of sound, 
the pressure and the mass densities are increasing in the 
interval (i?o,i?oo) : 

■j^X > for X = a 2 ,p and po- 

Proof of the Lemma 2. 

The estimate i?o *C Rao follows immediately from 
the boundary condition a 2 ^ 3> m/R^, and the assumed 
value of A. In the asymptotic end {RojRoa) the function 
—C(R) can be bounded from below as follows, 



C(R) 



-m(R) 

~~R 

-2m 



+ 2U 2 (R) 
AR 2 



4irR 2 p 



AR 2 



> 



(17) 



The second line follows from U 2 > and from the obser- 
vation that p < p — thus AnR 2 p < Ait fj^°° drrp — and 

4*f£~drrp<% 

Taking into account the last estimate, we can estimate 
from below the right hand side of Eq. ( 13 ) by 



D(R) 



1 



2m 



AR 
~3~ 



We obtain from Eqs ( 13 ) and ( 18 ) the inequality 

2m AR\ 



dR 



a 2 > (T-l-aL 



R 2 



3" 



(18) 



(19) 



It is clear that i) Ro is the null point of D(R); ii) the func- 
tion D(R) is strictly positive for R 6 (Ro,Roo)- There- 
fore a 2 (and consequently the remaining gas character- 
istics p and po ) does increase in this interval, provided 
that a solution does exist. 

The necessary condition for the existence of a solution 
is that a(Ro) should be strictly positive. Integration of 
Eq. (19), between i? an d -Rooi yields 



— a 
2m 
R 

2m 



2 (Ro) > (r - i - <4 

AR 2 2m AR 



We can replace by 
and get the following 



G 

3 



Ra 



G 



(20) 



using the definition of Ro, 



T-l-a 



, 2m o 
R — "oo 



r - 1 - at 



AR| 
18 



6x361 (itY 1 



(21) 



In the second inequality we dropped the positive term 
with m/Rrx, and used the definition od Ro- We already 
assumed that m/i?^ <C which ensures that the first 
term within second bracket is small and the right hand 



side of the last inequality in (21 1 is approximated by 
A " /r ' ~ 2 ). Therefore if the cos- 



E(R) 



(r-i 



mological constant is given by AR 2 ^ = 18a 



r-i- 



with 



a larger than 1, then E(R) must be negative and the 
square of the speed of sound at i?o must be nonposi- 
tive , which means that solutions are absent. On the 
other hand solutions can exist for 0.01 < a < 1; in this 
case C(R) < between (Rq,Roo), which implies that 



m x < for X 



a ,p, po and p. 



The two lemmas proven in this section have a trans- 
parent physical interpretation. Lemma 1 assumes that 
the dark energy density A /Air is smaller than one half of 
the averaged matter density 3m/47ri?^ . The presence of 
A is expected to have some quantitative impact onto ac- 
cretion mass rate M, but qualitative features of the flow 
are not influenced — in particular the gas characteristics 
a, |J7|,po an d p are all decreasing functions. 

Lemma 2 assumes the opposite — that the dark energy 
density A/47T is much larger than the average matter den- 
sity 3m/47ri?^ . Above a particular value — that depends 
on the boundary characteristics of the flow, its volume 
and the black hole mass — steadily accreting solutions 
do not exist. This is intuitively understandable — large 
A implies a large Hubble expansion velocity which can 
obstruct or even prohibit accretion, through the junction 
condition on the boundary of the vacuole. Therefore one 
can expect a significantly smaller mass accretion rate M 
and even its absence. 



QUALITATIVE ANALYSIS OF SOLUTIONS: 
SADS BLACK HOLES 



Lemma 3. Assume A < 0. Then: 
i) for R > R* the speed of sound, the pressure and 
e mass densities po and p are decreasing, j^X < for 
= a 2 ,p and po- 



< -6 



r-i- 
r-i- 



ii) Solutions are absent if AR 2 ^ 
Proof. 

Part i) The key observation is that the function C is 
strictly positive, if A < 0. Indeed, let Rb > R* be a small- 
est radius at which C(Rb) vanishes, and assume C < 
in the open interval (i?*,i?h) (notice that C(i?*) van- 
ishes). But Eq. (15 1 would imply in such a case that 



both terms with derivatives become negative in an sub- 
interval ,(R* , i?i ) of (i?*,i?fc), that is — see the remark 
below ( 10 ) — Oro 2 < 0. But this in turn would con- 
tradict Eq. (fl3l — if C < then dud 2 > outside the 



sonic sphere. Thus we arrive to the contradiction in this 
sub-interval (R*, i?i) and the argument of continuity and 
the definition of Rb imply the positivity of C in the open 
interval (R*, Rb)- 

In the next step we show that C cannot vanish outside 
the sonic sphere. Indeed, (13 1 and ( |15| imply that g§ 
is strictly positive wheneverC 1 = for R > R*; that 
in turn yields (by the argument of continuity) the strict 
positivity of C outside the sonic sphere. Eq. ( 13 ) allows 
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one to conclude that the speed of sound is monotonically 
decreasing, and that implies also decrease of the mass 
density po an d the pressure. 

Part ii) In order to prove the necessary existence con- 
dition, notice that Eq. ( 13 1 gives, after integration, 



|aL-« 2 CR)l > 



dr\ (r - 1 



-2U 2 



~Y 

(22) 

Notice that a 2 < T - 1; thus - a 2 (R)\ < T — 1 — 
aL, . The conservation of the mass accretion rate M and 



jfiPo < lead to the estimate 
in turn implies inequality 



d 

dR 



U 2 R 4 ) > 0, and this 



U 2 < Ul 



i? 4 



(23) 



Collecting this information, one finally obtains 



r-i-<&> (r-i- a 2 (i?)) - 



u 2 



AR 2 



G 



1 - 



If 
(24) 

AR 2 T— 1 — a 2 

Let R = R* and — < — r _ 1 _ ( ^2 ; then — taking into 
account that j=pr is negligibly small, and the same is true 

for the velocity term U 2 — we conclude that the right 
hand side of Eq. (22) exceeds its left hand side. This 



contradiction proves the necessary condition of Lemma 
3. 



VI. DESCRIPTION OF NUMERICS 

Our aim is to find transsonic accretion flows. We 
assume the total mass m — 1, T — 4/3 and Roo — 
10 6 . The boundary values at R^, are k — N — 

\A ~ IS ~ U ™ ~ t^D ^ ^ = 2 X 10 ~ 4 ' There 

three remaining parameters that are free: the cosmo- 
logical constant A, the mass accretion rate M and the 
asymptotic density pooo — after specifying them, equa- 
tions are integrated inward, starting from R^. These 
quantities should be such as to ensure the asymptotic 
conditions: \Uoo \ -C y / 2m/R oc> -C Ooo. 

A numerical run, with the above boundary data, has 
been regarded as succesful, if: 

i) the integration started from i?oo and reached the 
outermost apparent horizon, that the point Ri such that 
2m j^ R < ) + = \ t along the accretion branch (with the 
increasing — in the direction to the centre — infall veloc- 

ityjf/|); 

ii) there appears a sonic point at a radius i?* > Ri. 
More specifically, in order to find a single transsonic 

flow, we fix A and the mass accretion rate M, and treat 
the asymptotic value of the baryonic mass density as a 
free parameter. For a random choice of pooo there are 
three possibilities, one exceptional and two generic: 

a) an almost unlikely event, that the numerical run 
is succesful in the sense defined above and we find the 
sought accretion flow; 



b) the solution numerically fails — a singularity ap- 
pears during integration; 

c) there exists a numerical solution without a sonic 
point — a piece of the solution describes accelerating 
accreting flow, while the other corresponds to a deceler- 
ating accreting gas. 

In the case a) the job is done. If one of the remaining 
generic alternatives occurs — say b) , with the asymptotic 
density — we change pooo until we get the situation de- 
scribed in c) with the asymptotic density p c )ooo- Now we 
use the bisection method in order to find an intermediate 
value piooo between the two values Pb)ooo an d p c )ooo> f° r 
which the numerical run is succesful and the transsonic 
flow exists. 

The numerics itself is standard — we use commonly 
known algorithms of the 8th order. 

The calculations had been performed for a few dozens 
of values — negative and positive — of the cosmological 
constant A. In the sector of negative cosmological con- 
stant we take A e (-5 x 1CT 13 , -10~ 19 ) ; the absolute 
value over bound is larger by one order of magnitude 
than the value found in Lemma 3 as a necessary exis- 
tence condition. In the positive part of the spectrum we 
choose A e (l0~ 19 , 7.2 x 10~ 15 ); this is quite close to the 
value established as the necessary condition in Lemma 2. 
For each fixed value of A we find about 100 solutions 
corresponding to different values of M. 

The obtained information is summarized in figures 1 
and 2. The abscissa shows the ratio x — M g /m of the 

mass M g = 4tt J^°° drr 2 p of gas to the total mass m. The 
total mass has been normalized to 1. We obtain, although 
in a less transparent form, the same result as in ([5]) - 
the maximum of the mass accretion rate corresponds to 
i « 1/3. For a fixed value of M there exist two flows, one 
with heavier center and lighter gas, and the other with 
opposite characteristics. Just around the maximum of M 
there are blank points — we find that in this region it is 
difficult to find a numerical solution. We think that this 
a purely numerical artefact. There is a slight shift of this 
maximum towards smaller values of x with the increase 
of the absolute value of the cosmological constant. 

The only new relevant physically information is that 
concerning the mass accretion rate M. It appears that 
the presence of the cosmological constant manifests only 
when the dark energy density A/ (Aw) exceeds the aver- 
aged matter density 3m/(AwR^ >0 ) «4x 10~ 19 by at least 
three orders of magnitude. The mass accretion rate M 
is decreasing with the increase of the absolute value of 
the cosmological constant. It is clear that the rate of the 
falloff depends on the sign of A — it is faster for positive 
values — but the falloff itself occurs both for positive 
and negative cosmological constant. 



VII. COSMOLOGICAL IMPLICATIONS 

Results that are reported in the preceding section 
demonstrate that while the mass accretion rate depends 
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FIG. 1: The ordinate shows the mass accretion rate M and 
the abscissa shows 1 — x, where x is the relative mass of gas 
in the system.. The various lines correspond to Afi^/lO -3 = 
1, 5, 6, 7, 7.15 in the order from the top. 
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FIG. 2: The ordinate shows the mass accretion rate M and 
the abscissa shows 1 — a;, where x is the relative mass of gas 
in the system. The various lines correspond to Ai?j^/10~ 3 = 
— 1, —25, —50, —100, —250, —500 in the decreasing order from 
the top. 



on the cosmological constant, this effect becomes signif- 
icant only when the dark energy fraction Q\ is much 
larger than the material fraction Q m . On the other hand, 
this impact can be dramatic, a sevenfold increase of A — 
starting from its value corresponding to ilj^/rt m « 10 3 
— leads to the diminishing of M by seven orders in mag- 



nitude, as seen in Fig. 1. This is so for the accretion 
onto SdS black hole, but the phenomenon is strong also 
in the case of SAdS (Fig. 2) . When applying these re- 
sults to our Universe, we have to take into account, that 
3> in two epochs — during inflation and after 
10 11 years from Big Bang. 

Inflation strictly excludes steady accretion onto pri- 
mordial black holes, due to Lemmas 2 and 3, because 
during inflation era J7a exceeds £l m by something like 50 
orders of magnitude. This fact is concordant with the 
well known freezing of structures that are bigger than 
a particle horizon |18j . but the difference is that now it 
applies to small accretion systems that are located well 
within the horizon. 

In the present Universe the dark energy and material 
densities are roughly the same, which suggests that the 
effect of dark energy is negligible. Notice, however, that 
the tenfold ageing of the Universe would increase the ra- 
tio flj^/U m by a factor of 100. That means that in the 
Universe older than 10 11 years the steady accretion would 
become less efficient and at the time t — > 10 12 years its 
efficiency goes to zero. 

Penrose conjectured the so-called Weyl curvature hy- 
pothesis [inj. It asserts, in its informal version, that a 
Friedman-Lemaitre-Robertson- Walker spacetime, where 
the Ricci curvature is nonzero but the Weyl curvature 
vanishes, evolves towards a vacuum spacetime filled with 
black holes and gravitational radiation — with nonva- 
nishing Weyl curvature and negligible Ricci tensor. (A 
more formal scenario has been worked out by Tod |20) ) . 
This picture assumes that initially small inhomogeneities 
of FLRW universe accrete matter and transform them- 
selves into black holes, that gradually merge, leaving at 
the end a net of huge black holes and a lot of gravita- 
tional radiation. A cautious intepretation of our findings 
would be to say that the role of steady accretion in the re- 
alization of this scenario, in the presence of dark energy, 
is insignificant. It might well happen, however, that this 
property of dark energy of diminishing the steady spher- 
ical accretion signals a more general feature — that dark 
energy damps any accreting process. In such a case there 
arises another fundamental question: how dark energy 
impacts this scenario outlined in the Weyl curvature hy- 
pothesis? 
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